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ABSTRACT 

Dual AGNs are natural byproducts of hierarchical mergers of galaxies in the ACDM cosmogony. 
Recent observations have shown that only a small fraction (~ 0.1% — 1%) of AGNs at redshift z < 0.3 
are dual with kpc-scale separations, which is rather low compared to the high merger rate of galaxies. 
Here we construct a phenomcnological model to estimate the number density of dual AGNs and 
its evolution according to the observationally-estimated major merger rates of galaxies and various 
scaling relations on the properties of galaxies and their central massive black holes. We show that 
our model reproduces the observed frequency and separation distribution of dual AGNs provided that 
significant nuclear activities are triggered only in gas-rich progenitor galaxies with central massive 
black holes and only when the nuclei of these galaxies are roughly within the half-light radii of their 
companion galaxies. Under these constraints, the observed low dual AGN frequency is consistent 
with the relatively high merger rate of galaxies and supports the hypothesis that major mergers lead 
to AGN/QSO activities. We also predict that the number of kpc-scale dual AGNs decreases with 
increasing redshift and only about 0.02%-0.06% of AGNs are dual AGNs with double-peaked narrow 
line features at redshifts of z ~ 0.5 — 1.2. Future observations of high- redshift dual AGNs would 
provide a solid test for this prediction. 

Subject headings: black hole physics-galaxies: active-galaxies: interactions-galaxies: nuclei-galaxies: 
Seyferts-quasars: general 



1. INTRODUCTION 

Merger of galaxies is one of the major processes of 
galaxy formation in the hierarchical ACDM cosmogony. 
Galaxy mergers can naturally lead to the formation of 
paired and binary massive black holes (MBHs), since 
most galaxies, especially those with sp heroidal compo- 
nents, host MB Hs at their centers (e.g.. iBegelman et al.1 
198(153 [III). If both of the merging galaxies are gas 
rich, a large amount of gas can be channeled to the cen- 
tral region of each mergi ng galaxy, as in dicated by nu- 
merical simulations (e.g.. lHernquistl ll989i). A dual AGN 
(dAGN) could then emerge if the accretion onto both 
MBHs is triggered during the merging process. Thus the 
existence of dAGNs in the universe and their demography 
provide an important probe not only of the hierarchical 
galaxy formation models but also of the triggering mech- 
anisms of nuclear activities and the assembly history of 
MBHs. 

In the past decade, substantial progress has been 
made in searching for dAGNs and binary MBHs (BBHs) 
through various possib le signatures, such as double- 
peaked broad lines fe.g.. lBoroson fc Lauerl I2009TV double- 
peaked narrow lines (e.g~T|Zhou et al.ll2004t"lGerke et al.1 
2007t IComerford et al.l l2009al|bl: IXu fc Komossal 1200 ' 
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Wang et a 



1 120091 : iSmith et all 1201(1 iLiu et al.l l20"Toal[ 



Shen et all l2010t iFu et al.l 120101 : IRosario et al.1 l201li 

and various other methods (e.g. , iKomossa et al.l 7200 
iValtonenl 120081 : [Green et aT1l2010f ). Following a system- 
atic scrutiny of the NIR images and optical slit spectra 
of a sample of double-peaked narrow line type 2 AGNs, 
it has recently been reported that roughly 0.5%-2.5% of 
the z < 0.3 type 2 AGNs are kpc-scale dAGNs after tak- 



ing in to account the sele ction co mpleteness (ILiu et al.l 
I2010allbl : IShen et alj [2010L see also IRosario et al.l 120 111) . 

This frequency is surprisingly low even if we were to as- 
sume that all these observations are indeed of dAGNs and 
not as has been otherwise suggested of superposed sin- 
gle AGNs, bipolar jets and accretion disks. The observed 
frequency of dAGN is indeed more than an order of mag- 
nitude lower than that expected from the observed major 
merger rate of galaxies if one makes the as sumption that 
each major m erger would yield a dAGN (jRosario et al.1 
120111 see also IShen et "all 1201 Of) . The lo w frequency of 
dAG Ns, unlikely due to selection effect (jRosario et al.l 
1201 lh . seems generally at odds with the commonly ac- 
cepted scenario that AGN/QSO activities are triggered 
by major mergers of galaxies. 

In this paper, we construct a phenomcnological model 
to address the observed frequency of kpc-scale dAGNs, 
by taking into account the following factors: (1) the 
merger rate of galaxies; (2) the types of galaxy mergers 
that can trigger the nuclear activities of both progen- 
itors; (3) when and where significant nuclear activities 
are triggered (by "significant" nuclear activities we mean 
that the Eddington ratios of the nuclear luminosities are 
close to 1 (e.g., 0.1-1), rather than 10~ 3 or less); and (4) 
how long the nuclear activities of both progenitors can 
last before the merger of the two MBHs. After the merger 
of two MBHs, the nucleus of the merged galaxy is likely 
to still be active, but it only appears as a single AGN. 
In Section [2l we show how these factors are adopted in 
our model of the dAGN distribution. Under the assump- 
tion that significant nuclear activities can be triggered 
only in gas-rich progenitor galaxies with central MBHs 
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when their companion galaxies are sufficiently close to 
them, the dAGN frequency and their separation distri- 
bution are obtained in Section [3] By comparison with 
observations, we show that the observational frequency 
of dAGNs is consistent with current observational con- 
straints on the merger rates of galaxies and the scenario 
that major mergers of galaxies lead to significant nuclear 
activities. In Section [3l we also present our model pre- 
dictions for the dAGN distribution at higher redshifts. 
We further discuss the uncertainties in our estimates of 
the dAGN frequency due to different processes, which are 
not included in our model, such as mergers of progenitor 
galaxies whose nuclei have already been activated in pre- 
vious mergers and also tidally-induced nuclear activities 
in two galaxies which are still far away from each other. 
Conclusions are given in Section 2] 

2. A PHENOMENOLOGICAL MODEL OF DUAL AGNS 

In this section, we introduce a phenomenological model 
to estimate the frequency of dAGNs, by taking into ac- 
count the factors mentioned in Section [TJ The crucial 
elements of this model are described as follows. 

2.1. Merger rates of galaxies 

For galaxies with stellar mass > M* , we describe their 
merger rate at a given time by the number fraction of 
those galaxies for which a galaxy merger completed per 
unit time, and denote the merger rate by 



slope f3(M t ) is - 0.15 - 0.20. The dependence of the 
merger rates on mas s ratio x can be approximated by 



1l{> M*,x,z) 
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dn mrg (> M*,x,z) 



«tot(> M*,z) 



dt. r 



(1) 



where t z = J°° \jp-\dz' is the cosmic time at redshift z, 
x(< 1) represents the mass ratio of two merging galaxies, 
fitot(> M*, z) is the comoving number density of galaxies 
with stellar mass > M* at redshift z, and (dn mrg / dt z )dt z 
gives the comoving number density of those galaxies that 
are products of mergers of two progenitor galaxies with 
mass ratio > x and with their mergers being completed 
over a cosmic time from t z to t z + dt z . 

In the past several years, tremendous efforts 
have been made in estimating the merger rate 
of ga laxies either through c lose pairs of galaxies 



e.g., iKartaltepe et ail 120071: iLin et all 
de Propris et al.ll2005t iBell et all l2006allbl: 
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20M. [20081: 



Bundv et al.1 



20091 lLopez-Saniuan et all l200aTde Ravel et alj 120091: 
Patton et al.l |2011|) or through morphological distur- 
bances of galaxies found in various deep surveys (e.g. , 
Cassata et al.ll2005t iBundv et al.ll200l IWolf et al.ll2005l: 
Bridge et al.ll2007ll2010HLotz et al. 2008; Conselice et"ai1 



20081120091 ). The dependence of the merger rate TZ(> 



M*, x, z) on mass ratio x can be absorbed in a function 
f(x), which describes the fraction of mergers with mass 
ratio larger than x, i.e., 

1Z(> M», x, z) = f(x)R(> M*, z). (2) 

For major mergers (usually defined by a; > 1/3), observa- 
tional estimates can be roughly fitted by the simple for- 
mula, R(> M„z) = A(M*) major (l + z)P( M *^i°*, where 
A(M») maj o r « 0.2[1 + (Ah/AI ) a - 5 } Gyr" 1 , /3(M*) major a 
1.65 - 0.151og(MJ Af n ), and M = 2 x 1O 1O M 
([Hopkins et al.ll2010t ). The normalization of the merger 
rate A(M*) is uncertain by a factor of about 2 due to 
systematic errors and the uncertainty in the evolution 
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on mass rati o x can be ap proximat ed by 
(1 - x) ([Stewart et al.l 120091: Hopkins et all 

[2010h . 

The stellar mass function of galaxies, involved in equa- 
tion ([T]), has been estimated over a large redshift range 
from various galaxy redshift surveys. In our Monte- 
Carlo calculations below, we adopt the latest stellar mass 
function obtained from the S-COSMOS survey (Spitzer- 
Cosmic Evolution Survey) for galaxies at redshift ranges 
0.2 - 0.4, 4 - 6, 0.6 - 0.8, 0.8 - 1.0, and 1.0 - 1.2 
([Ilbert et al.1 12010| ) , and that obtained from the Sloan 
Digital Sky Sur vey (SPSS) for galaxi es in the local uni- 
verse (z ~ 0.1) ([Bernardi et al.ll2010D . respectively. Ac- 
cording to equation ([!} and the stellar mass functions, 
our simulations generate a large number of merging pairs 
of galaxies over redshift to 1.2. The masses of the 
two progenitors are assigned by M*^ = max(^, 

and M* t 2 = mm (f^, iq^), respectively, where x is 
randomly selected according to its distribution function 
fix)- 1 

2.2. Triggering of nuclear activity and mergers of 
galaxy pairs with various morphologies 

During the merging process of two galaxies, whether 
significant nuclear activities at their centers can be trig- 
gered depend on the two factors: (1) whether MBHs are 
initially hosted in both nuclei or whether the MBHs are 
massive enough; and (2) whether sufficient gaseous ma- 
terials can be quickly delivered into the vicinity of the 
MBHs. Those two factors are closely related to the pro- 
genitor morphologies of the merging pairs. 

2.2.1. The mass of the initial MBHs 

Observations have shown that the masses of MBHs, 
M,, in the centers of nearby galaxies are tightly corre- 
lated with the stellar ma ss of the spheroidal components 
of the galaxies A/^buigc 

([Haring fc Rixll20"ol 
log M. = 8.28 + 1.12(IogM.,bnige - 11), (3) 

with an intrinsic scatter of 0.3 dex in logM,. Subse- 
quently, given the stellar mass of the spheroidal compo- 
nent M^buigc.i (i = 1, 2) of each progenitor of a merging 
pair, we use equation ([3]) to estimate the central BH mass 
M, t i (i = 1,2) in the progenitor at redshift z, by assum- 
ing that the scatter of the correlation follows a normal 
distribution and adopt ing a small evolution correction of 
M. oc (1 + zf m (e.g- lMerloni et al.ll2010h . 

Note that the following elements are involved to esti- 
mate the distribution of M^buigev 

(a) The bulge to total stellar mass ratio (B/T): given 
the stellar mass M* of a galaxy, the mass M^buige can 
be estimated according to the morphology of the galaxy, 
as the B /T ratios a re different for galax ies with dif- 
ferent morphologies ([Weinziri et al.1 l2009f ). Based on 
the detailed analysis of a sample of nearby galaxies, it 

1 Note that a substantial number of stars may be formed in 
the merging systems during the merging process. We have checked 
that if a fraction (e.g., 10% — 50%) of the total stellar mass M, 
is removed in order to get the real total stellar mass of the two 
progenitor galaxies before doing the splitting procedure above, a 
slightly lower estimate of the dAGN frequency may be obtained, 
but which does not affect our conclusions in Section [4] 
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has been found that B/T = 0.22 and 0.05 with vari- 
ances of 0.05 and 0.0 2 for Sa-Sb and Sc-Sd, respectively 
(|Weinziri et al.ll2009t ). We assume that B/T =1 for el- 
liptical and SO galaxies and B/T = for irregular galax- 
ies, respectively. 

(b) Distribution of morphological combinations of 
merging galaxy pairs: it is possible that the two progeni- 
tors of a merging galaxy pair have different morphologies. 
Currently detailed statistics on the morphological combi- 
nations of the pairs are not available, although there are 
quite a number of constraints on the merger r ates of red 
galaxies, b lue galaxies and mixed ones (e.g., iLin et ail 
12004 120081 ). Here we simply assume that the morpho- 
logical types of the two progenitors are independent of 
each other, and consequently the fraction of the merger 
rates of the pairs with any specific morphological com- 
bination to the total merger rates only depends on the 
fraction of these types to the total number of galaxies. 
This assumption appears consistent with the current ob- 
servations that the fractions of different morphological 
types for galaxies in pairs are similar to those for field 
galaxies (|de Propris et al.ll2005D . 

(c) Stellar mass functions of galaxies with different 
morphological types: the galaxy luminosity function 
(LF) for four spectral types of galaxies over z ~ 0.2 — 1.2 
was previously estimated (see Table 3 of IZucca et al.1 
2006). These four spectral types roughly correspond 
to the morphological types E/S0, Sa-Sb, Sc-Sd and Irr, 
respectively. For each type of galaxies, the mass-to- 
light ratio can be estimated through their average col- 
ors (e.g.. B-V. seelBell et al.ll2003t feernardi et al l 20101: 
iFukugita et al lll995l) . and their LFs (IZucca et alj|2006ft 
can thus be converted to the stellar mass functions. Ac- 
cording to these mass functions, the fraction of each type 
of galaxies to the total galaxies can be obtained at any 
given M* over z ~ 0.2 — 1.2. For galaxies at redshifts 
z < 0.2, we adopt the stellar mass functions for diffe rent 
morphological types given bv iBernardi et al.l (|201d . sec 
their Table B2) and estimate the fraction of each type to 
the total at any given M*. In our calculations below, for 
each merger generated as that described in Section |2~T1 
we randomly set the morphological type to each progen- 
itor according to the fractions of each type of galaxies to 
the total at M*^ and M* j2 , respectively. 

For those irregular or purely disk-like galaxi es, we set 
the initial mass of the ir MBHs to be 10 5 M o |Hc][2008; 
iKormendv et al.ll2011l ) or zero. Here, these two different 
initial values do not lead to significantly different results. 

2.2.2. Triggering of the nuclear activity 

During the merger of late-type gas-rich galaxies (Sa- 
Sb, Sc-Sd, or Irr) with other galaxies, significant nu- 
clear activities may be triggered rapidly. The reason 
lies in the fact that late-type galaxies contain signifi- 
cant amount of gaseous materials which could lose angu- 
lar momentum and sink under dynamical friction during 
the galaxy merger. Early-type gas-poor (red/elliptical) 
galaxies may be able to capture some gas from gas-rich 
encounters during the merging process, but the time for 
the gas to reach the center is likely to be comparable 
to the merging timescale (see Section |2.3[) . Therefore, 
the nuclear activities in these galaxies may start only 
after the mergers have almost completed. In addition, 
the host galaxies of all the confirmed dAGNs selected 



through double-pe aked narrow lines contain stellar disks 
(|Shen et al.l 12010ft . which would suggest that they are 
merging remnants of late-type galaxies. 

Based on the above arguments, we assume that the 
nuclear activity is triggered in each gas-rich component 
with a central MBH once its separation with its compan- 
ion becomes smaller than a threshold D c . On the other 
hand, we assume that the nuclear activity cannot be trig- 
gered in any gas-poor component of a merging pair prior 
to the completion of the merger (see also discussions in 
Section [3]). The physical size of a galaxy is characterized 
by its half-light (or effective) radius r^. Once the sepa- 
ration of the component i (1 or 2) of a pair to the other 
component j (2 or 1) is smaller than the half-light ra- 
dius of component j (i.e., rh.j), the center of component 
i may be significantly perturbed. Here we assume the 
threshold D c for component i to be D c i = Kr^ j, where 
K is a fudge factor of order unity. For major mergers, 
nuclear activities can be triggered in both components 
only if both components are gas-rich and their separa- 
tion D < D c = Xmin(rh,i, ^,2)7 where 7"h,i and fh,2 
are the half-light radii of the primary and the secondary 
galaxy of the pair, respectively. 

For disk galaxies, we adopt the following empirical re- 
lation to estimate the disk half-light radius from the stel- 
lar mass, which is obtained from a SPSS sample at low 
redshift z ~ 0.1 (jDutton et al.ll2011l ) 
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where a 
and Mq 
lation at 
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1.8, r = 5.2 kpc 
The scatter of this re- 
to follow the log- 



= 0.18, (3 = 0.52, 
= 2.75 x 1O 1O M . 
a given M* is assumed 
normal distribution with a standard deviation of a 
S2 + (si-s 2 )/[l+(M*/Mo)"], where Sl = 0.47, s 2 = 0.27, 
Mq = 2 x 1O 1O M and 77 = 2.2. For disk galaxies at 
higher redshifts, for example z ~ 0.2 — 1.2, we add an 
evolution correction to the zero point of the disk-size — 
stellar mass relation by A log ro — 0.018 — 0.44 log(l + z) 
(jPutton et al.ll2011h . 

For early- type galaxies, the half-light size versus stellar 
mass relation can also be fitted by equation (j4]), but with 
a = 0.03, /3 = 0.64, 7 = 1-3, r n = 1.4 kpc, and M = 
1.2 x 10 10 M o (jPutton et al.ll2010D . This double power- 
law form is consistent w ith the single pow e r-law form 
at high masses g iven by Ivan der Wei et all (|2008T ) and 
IShen et all (|2003h . The intrin sic scatter of this relatio n 
at high masses is ~ 0.14 dex (Ivan der Wei et afll200ll . 
We assume the scatter is also ~ 0.14 dex at low masses 
(<~ 10 9 — 10 10 M Q ). For early-type galaxies at higher 
redshifts, for example z ~ 0.2 — 1.2, we add an evolution 
correction to the zero point of this relatio n as A log ro = 
-0.981og(l + z) (|van der Wei et al.l 120081 ). 

2.3. Evolution of the separation of two merging galaxies 

The merging timescale of paired galaxies is crucial in 
determining the time period during which the merging 
system would appear as dAGNs. Based on virtual galaxy 
catalogues obtained from the Millennium Simulation, it 
has been found that the average merging time of paired 
galaxies can be fitted by 

\ -0.3 



(T mrg )(r p ) = 2.2Gyr 



50 kpc/ V4 x lO 10 ^- 1 .^ 
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(5) 



where r p is the projecte d separation of the pairs 
(|Kitzbichler fc White! [20081) . This result is about 15% - 
30% larger than that obtaine d from hig h er res olution 
simulations of galaxy mergers lLotz et all (|2010[ ). Cor- 
recting this difference and assuming that the orientations 
of galaxy pairs are isotropically distributed, the merging 
time as a function of three-dimensional separation D is 
smaller than that given by equation ([5]) for a projected 
separation r p = D by about a factor of 2, i.e., 

(r mrg ) (D) ~ 0.5 (T mrg ) (r p = D). (6) 

The period that a merging system may appear as a 
double-peaked narrow line dAGN varies from St = to 
5t = T mrg (D c ) — r mrg (Z?L), where £>l is the lower limit 
of the separation that the narrow line regions associated 
with the two nuclei begin to overlap and the line profile 
becomes more complicated, and 5t is the time elapsed 
since the activities of both the nuclei were triggered. 
Here we set = 0.5 kpc since the typical size of nar- 
row line regions of AGNs is about a few tens to a few 
hundreds of parsecs (|Petersonlll99l . 

2.4. Luminosity evolution of nuclear activities 

Once the nuclear activity is triggered in the center of 
a merging component, the material available for accre- 
tion initially is expected to be abundant. The situa- 
tion of these merging systems may be similar to that 
of high luminosity QSOs, which also have plentiful gas 
supplies and accrete at relatively high rates. The Ed- 
dington accretion rate is MEdd — 0.22A/ Q yr _1 (e _1 — 
1)(M./1O 8 M ), where e ~ 0.1 - 0.16 is the mass-to- 
light con version efficiency (e.g., lYu fc Tremamel 120021 : 
lYu fc Ldl2008h . Here we assume that the accretion rate 
is AMEdd for the MBH in each component of a merg- 
ing pair after its nuclear activity has been triggered, 
and we set A = 0.25 (e .g., as shown in iKollmeier et all 
I2006t IShen et all l2008ft unless otherwise stated. Thus 
the MBH should grow exponentially with the elapsed 
time as exp(A<5i/rsai), where the Salpeter timescale rs a i 
is about 4.5 x 10 8 (Y^)yr. We assume that the accre- 
tion switches off if the total mass of the two MBHs is 
larger than the expectation from equation ([3]) for the final 
merged galaxy, taking into account the scatter in the re- 
lation. 2 The bolometric luminosity is Lboi = AiEdd(A/.), 
where L E dd ^ 1-3 x 10 38 erg s -1 (M./M ). For the 
purpose of this paper, we need calculate the [OIII] lu- 
minosity, which can be converted from the bolomet- 
ric luminosity by £[ in] ~ iboi/3500 with a scatter of 
~ 0.38 dex (assumed to have a log- normal distribution) 
(jHeckman et al.ll2004l) . 

2.5. Velocity separation 

In recent systematic searches, double-peaked narrow 
lines have been adopted as the indicator to search for 
dAG Ns from parent SDSS AGN samples (jShen et al.l 
20 10:. In this search the dAGNs with velocity separa- 
tion v sep < 150 kms -1 or v sop < max(er Ci i, a c ^) cannot 

2 Even if these BHs can still accrete some material via lower Ed- 
dington ratios, e.g., ~ 0.001, our simulation results are not affected 
much. 



be detected. The velocity dispersions of the primary and 
the secondary progenitor galaxies (or their bulges) of the 
pair, ct c ,i and <7 Ci 2 respectively, characterize the widths 
of the two narrow line components. In order to extract 
the dAGN frequency from these searches, it is necessary 
to estimate v sep , a c .i and cr Gj 2 of the two merging com- 
ponents. 

As both components of a pair move in their common 
dark matter halo with a velocity of the order of the cir- 
cular velocity v c i rc of the halo (i.e., the merged galaxy), 
we first assume that the relative velocity of the two com- 
ponents w ro i ~ 2v c i rc and later discuss the consequences 
of this assumption in Sectional It is plausible to assume 
that the orbital orientation of the two merging compo- 
nents and their relative position on the orbital plane are 
isotropically distributed. Subsequently the velocity sepa- 
ration of the two components of the double-peaked lines, 
i.e., the relative velocity component projected onto the 
line of sight, as well as their projected separation perpen- 
dicular to the line of sight, can be taken into account for 
each pair through Monte-Carlo simulations. For early- 
type galaxies, t> c j rc c an be estimated from u c i rc ~ 1.54er c 
(jDutton et al.l l2010h . The velocity dispersion a c of the 
merger result which is supposed to be an elliptical, can 
be estimated from stellar mass M* by the following em- 
pirical relation 



log 



kms 



= 2.23 + 0.37 log 



-0.19 log 



1O 1O9 M 
A/* 



10 ia9 M, 



(7) 



for nearby SDSS gala xies which has an int rinsic scatter 
of roughly 0.071 dex (jGallazzi et al.ll2006f ). For galax- 
ies at higher rcdshifts, we add an evolution correction 
to the above relation by in terpolating its mea sured de- 
viations to higher rcdshifts (jDutton et al.ll201ll see their 
Table 3). For a late- type galaxy, the velocity dispersion 
of its bulge (if any) can also be estimated from equation 
(0 after replacing M* there by M^buige, as bulges ap- 
pear to fol low the same Faber - Jackson relation as faint 
ellipticals (jBalcells et al.ll2007[) . 

2.6. Monte-Carlo simulations of dAGNs 

We summarize the procedures to generate dAGNs in 
our Monte-Carlo simulations as follows. 

• We generate 10 7 major mergers of galaxies which 
are completed during the cosmic time t Zk to t Zk + 
5t Zk , where z k = 0.1,0.3,0.5,0.7,0.9, and 1.1 for 
k = 1,2,..., and 6, respectively, and St Zk = 10 9 yr. 
The detailed completion time t z of a major merger 
is selected according to the major merger rates and 
the stellar mass function given in Section [2~TI For 
each merger, the total stellar mass M* (1O 9 M — 
10 12 Mq) of the system is split into the stellar 
masses of its two progenitor galaxies M*^ and 

• For each merging system, we assign morphological 
types to the two progenitor galaxies and then es- 
timate their bulge masses and the masses of their 
central MBHs, as described in Section [2. 2. II 
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Fig. 1. — The LFs of AGNs and the bright components of dAGNs. 
The red and blue dashed lines represent the observed [OIII] LF of 
all (includin g type 1 and type 2) AG Ns at redshifts z = 0.1 and 0.3, 
respectively (Bongiorno ct al. 2010). The dotted lines are reference 
lines. The red dotted line represents 1% of the observed AGN 
[OIII] LF at z = 0.1, and the blue dotted line represents 0.1% of 
the observed AGN LF at z = 0.3. Each nucleus of a dAGN has its 
own luminosity, and thus the luminosity of a dAGN is composed 
of two components. The red and blue solid lines represent our 
model-predicted LF of the relatively bright components of dAGNs 
at 2 = 0.1 and 0.3, respectiv ely. Our simulated dAGNs are selected 
by the criteria discussed in 5 12.51 The typical error to the estimates 
of the dAGN LF shown in this Figure is 0.5 — 0.6 dex. 

• For each progenitor galaxy of the system, the 
half-light radius is estimated as described in Sec- 
tion 12 2.2[ and their separation threshold for trig- 
gering nuclear activities is correspondingly set to 
each component, i.e., D C i = Kr^j (i,j = 1,2 and 
i 7^ j), where K is assumed to be the same for the 
two components of the major merger and is cali- 
brated by the separation distribution of observed 
dAGNs (see Section HJ}. 

• We start from the completion time of each major 
merger and trace back the evolution of the sepa- 
ration of its progenitor galaxies from zero to the 
values at earlier time, as described in Section 12.31 
(Note that our analyses for kpc-scale dAGNs here 
are not affected by our extrapolation of equation 
© to pc scales and by our simplification of de- 
tailed evolution timescales of massive bound binary 
black holes on these small scales.) Assuming that 
all merging systems which complete their mergers 
within the cosmic time t Zk to t Zk + 8t Zk are located 
at the distance of redshift z^, then these systems 
should appear in observations at a time of At = 
tz~t Zk before the completion of the mergers, where 
t z is the cosmic time at which each major merger is 
completed. For any merger, if At < (r mrg ) (-D c ,i), 
the nuclear activity in component i has been trig- 
gered, and the time elapse since the triggering of 
the activity is given by (Tmrg) (-Dc.i) — At. 

• The luminosity is assigned to each active compo- 
nent as described in Section [2~4l 

• If the nuclear activities in both components are 
triggered and the system satisfies the velocity 



threshold set in Section 12.51 it then appears as 
a dAG N similar to those observed by IShen et all 
d20ioh . 

In order to compare with the results given by system- 
atic surveys by using double-peaked narrow lines (e.g., 
IShen et al]|2010t iRosario et al.ll201l"l) . we only count the 
number of those simulated dAGNs with two components 
having a separation within the range from 0.5 kpc to 
lOkpc and comparable luminosities, i.e., the luminosity 
ratio of the two components is within a factor of 4 here. 
And we then calculate their luminosity function and pro- 
jected separation distribution, etc., as illustrated in the 
following section. 

3. MODEL RESULTS AND DISCUSSIONS 

Figure [T] shows the [OIII] LF of the bright components 
of the simulated dAGNs that could be selected through 
double-peaked narrow lines (solid and long-dashed lines) 
at redshifts z = 0.1 and 0.3, respectively. These sim- 
ulated dA GNs are se l ected through similar thresholds 
as that in Sh cn et ail (12010D for double-peaked narrow 
line dAGNs (see Section I2.5[) . And the threshold on 
the separation for the nuclear activity to be triggered, 
D C i = A>h.j, is set to have K = 1.25 (which is the 
reference value, see also Figure [2]), where £ = 1,2 repre- 
sent the two components of a merging pair and i ^ j. 
Note that the LF of the faint components of dAGNs 
is only slightly smaller than that for the bright compo- 
nents. The observed [OIII] LFs for all AGNs, including 
the [OIII] LF s of the type 1 AGNs and type 2 AGNs 
adopted from iBongiorno et al.l (|2010f ) , are also shown in 
Figure [T] Apparently, the LF of the bright components 
of dAGNs is smaller than that for all AGNs by a factor 
of about 100 to 1000. Note that the bright component 
of an observed dAGN may be h osted in the small pro- 
genitor galaxy (|Shen et al.ll2010t ): and our model shows 
that about 30% - 40% of simulated dAGNs have their 
bright components being hosted in the progenitors with 
relatively low total stellar masses, consistent with cur- 
rent observations. Here we do not try to simultaneously 
fit the total AGN LF, as both the detailed luminosity 
evolution of nuclear activities in the post-merger stages 
and other mechanisms leading to nuclear activities would 
be important in such a fit but beyond the scope of this 
paper. 

The accuracy of our estimates of the dAGN LFs de- 
pends on the uncertainties in (1) the merger rate TZ, (2) 
the merger timescale (T mrg ) (D), (3) the setting of the 
separation threshold D c , (4) the setting of the accretion 
rate A, and (5) the approximation of the relative velocity. 
We discuss each of the points below. 

• The systematic error of the merger rates, estimated 
from observations, could affect our estimate by 
a factor of 2. This error partly r eflects the un- 
certa inty in the merger timescale (|Hopkins et al.l 
I2010Q. The minor mergers of galaxies, not included 
in our calculation, might also contribute to dAGNs. 
However, the rate of minor mergers with mass ratio 
x £ (1/10,1/3) is almost the same as the rate of 
major mergers (jHopkins et alJl2010T >. The contri- 
bution by these minor mergers should not be larger 
than a few tenth of the contribution by the major 
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merger of galaxies (with x G (1/3,1)), as the two 
merger components differ substantially in stellar 
mass and are less likely to contribute to AGN pairs 
with two comparably bright components. There- 
fore, the total error introduced by the uncertainties 



in 1Z and (r n 



(D) is roughly a factor of 2. 



• The uncertainty in the settings of the separation 
threshold D c may be calibrated by the observa- 
tional separation distribution of dAGNs. Figure [2] 
shows the cumulative distribution of the number 
of dAGNs versus the projected separation r p of its 
two components for both observed dAGNs and sim- 
ulated dAGNs. As seen from Figure Eta), the sim- 
ulated distribution matches well the observational 
one if D Cji is set to ~ 1.25rh., . 3 After adding 
the dA GN candidate found by IComerford et al.l 
(|2009bl ) and the two dAGNs — major merg er sys- 
tems recently found bv lRosario et al.l (|201 lh to the 
dAGN sample, the dAGN r p distribution appears 
to be still consistent with D c ,i ~ (1 — 1.5)rhj. 
We choose K = 1.25 as the reference value in 
this paper. A much larger or smaller D c (e.g, 
K = 2 or 0.5) appears not to be consistent with the 
observationally-deduced distribution of projected 
separations (see also the discussion on tidal inter- 
actions below). According to equation J5]), the un- 
certainties in the settings of D c may lead to an 
error of about 30% in the simulated LFs of the 
bright (or relatively faint) components of dAGNs. 
Figure &b) shows the evolution of the dAGN sep- 
aration distribution. The dAGNs at higher red- 
shifts appear to have smaller separations, which is 
mainly due to the size evolution of galaxies set in 
the model. 



The accretion rates of dAGNs are set to have 
A = 0.25 times the Eddington rate, similar to 
those in observed bright QSOs (see Section f2.4j) . 
If A is set to be 0.1, the LFs of the simulated 
dAGNs move downward by a factor of 2 — 3 at 
£[Oiii] ~ 10 7 - 10 8 L Q ; and if A = 1, the LFs of the 
simulated dAGNs increase by more than a factor 
of 10 at I/[oiii] ^ 10 8 ' 5 Lq but stays the same and 
then decreases slightly as £[oni] decreases to less 
than 10 7 - 5 L G . 

The relative velocity v Te i of the two components of 
any merging pair is set to be twice the circular ve- 
locity of their common dark matter halo. This may 
be an overestimate of v ro i as the primary compo- 
nent may move at a lower speed than the circular 
velocity. However, it is hard to address the uncer- 
tainty in the dAGN LF caused by this assumption. 
In order to get a sense of this uncertainty, here we 



simply assume a case of v r 



l.5v c 



or v c 



and 



find that it could lead to a decrease of the dAGN 
LFs by about a factor of 1.4 or 2. 



3 W e caution here that the total number of dAGNs in lShen et al.l 
Il2010r ) sample is only 5, and the separation distribution of these 
dAGNs shown in Figuref2]suffcrs from small number statistics. Fu- 
ture observations which promise to find more dAGNs would provide 
firmer constraints on D c and would yield better statistics. 




r p (kpc) 

Fig. 2. — The cumulative distribution of projected separation 
of dAGN s. The histogram shows the observational results ob- 
tained by IShen et al.1 (2010). The simulated dAGN sample is se- 
lected by the luminosity of the bright component with i[oin] > 
lO 7,5 Z/0, as is set for the lowest luminosity threshold of the ob- 
served dAGNs. The simulated separation distribution is normal- 
ized to r p = 1.5 k pc, as is done for the observational distribution 
(Shcn et al. 2010). Panel (a) shows the results of our simulations 
at redshifts z = 0.1 (solid lines) and 0.3 (dotted lines), respec- 
tively. The separation threshold for triggering nuclear activities, 
i.e., D c i = Kr^ j (i = 1,2 and i ^ j for the two merging com- 
ponents) is set to K = 1 (red lines), K = 1.25 (blue lines) and 
K = 1.5 (green lines). Panel (b) shows the evolution of the distri- 
bution of the projected separations for the case of K = 1.25. The 
red solid, blue dotted, magenta dashed, cyan long-dashed, black 
dot-short-dashed, and green dot-long-dashed lines represent the 
results of our simulations at redshifts z = 0.1, 0.3, 0.5, 0.7, 0.9, 
and 1.1, respectively. 
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Fig. 3. — The expected dAGN frequency at different redshifts, 
i.e., the ratio of the number density of simulated dAGNs to the 
number density of AGNs derived from the observational AGN LF. 
Open circles and squares represent the cases for the bright com- 
ponents of dAGNs and AGNs with luminosities £[om] > 10 65 Lq 
and 10 7 ' 5 Lq, respectively. The dAGNs are selected through the 
double-peaked narrow lines in our Monte Carlo simulations de- 
scribed in Section \3\ In this figure, the separation threshold for 
triggering nuclear activities is set to have K = 1.25. The typi- 
cal error on the estimated dAGN frequency is 0.5 — 0.6 dex (see 
discussions in Section [3} . Th e total [OIII1 LFs of AG Ns at differ- 
ent redshift are adopted from Bongiorno et al. (2010). The dotted 
and solid crosses are the current observational constraints obtained 
with and without completeness correction, respectively. See details 
in Section [3] 
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Combining all these uncertainties together, the total 
systematic uncertainty in the LFs of the bright compo- 
nents of dAGNs is likely to be a factor of 3 — 4 over the 
luminosity range 10 6 — 10 s 5 Lq. 

Figure [3] shows the evolution of the dAGN frequency, 
i.e., the ratio of the number density of simulated dAGNs 
selected through the double-peaked narrow line features 
to the number density of AGNs derived from the obser- 
vational total AGN LF (|Bongiorno et al.ll2010h . Here we 
do not distinguish between type 1 and type 2 AGNs, 
but include both. This ratio is equal to the ratio of 
the cumulative LF of dAGNs to that of AGNs, and de- 
pends weakly on the lower limit of the AGN luminosi- 
ties. These estimates could be off by a factor of 3 — 4 
as previously discussed. At redshift z ~ 0.1 — 0.3, the 
estimated dAGN frequency is ~ 0.1% — 1%. (This value 
could be higher than the estimates from observed dAGNs 
selected through the double-peaked narrow lines, consid- 
ering that the signal-to-noise of some dAGN spectra are 
so low that the dAGNs arc not detected observation- 
ally.) There are quite a number of dAGNs which cannot 
be detected by using double-peaked narrow lines, due to 
the projection effect and the velocity cutoff. Including 
the undetected dAGNs, our model suggests that the real 
frequency of dAGNs is ^0.2% — 2%, larger than that of 
those dAGNs selected through the double-peaked narrow 
line features by a factor of ~ 2. These values are fully 
consistent with t he observational estimates obtained by 
IShen et al.1 (|201Cf ). i.e., - 0.1%-0.5% without complete- 
ness correction or ~ 0.5% — 2.5% after the correction of 
in-completeness due to the projection effect, velocity cut- 
off and low signal-to-no ise spectra of some dAGNs (see 
also lRosario et aill2011l in which the dAGN frequency is 
~ 0.4% without the completeness correction). 

The reasons that a lower frequency of dAGNs is ob- 
tained from our m odel (compared to the simple esti- 
mates; e.g., 15% bv lRosario et al-lfeOlll ) are three folds: 
(1) dAGNs can be produced only if both of the two pro- 
genitor galaxies are gas rich so that sufficient accretion 
materials could be provided; (2) MBHs in those gas- 
rich progenitor galaxies may be initially substantially 
smaller than their final masses, and thus the AGN phe- 
nomena triggered before MBH mergers may be substan- 
tially less luminous compared with those triggered after 
MBH mergers; and (3) in a significant fraction of major 
mergers, the masses of the MBHs hosted in the two pro- 
genitor galaxies may differ by a large factor and thus the 
induced dAGNs may have too large luminosity difference 
to be included. 

As shown in Figure El the modeled frequency of kpc- 
scale dAGNs that can be detected through the double- 
peaked narrow line features declines to even lower values 
on the order of several ten-thousandths, e.g., 0.02% — 
0.06% at redshift z ~ 0.5 — 1.2. And the frequency of 
all kpc-scale dAGNs declines to 0.04% - 0.1% at redshift 
z ~ 0.5 — 1.2. This decline is mainly due to the redshift 
evolution of the galaxy morphology distribution (e.g., 
IZucca et al.l 120061 ). Correspondingly, the major merg- 
ers at such redshifts are dominated by very late-type 
galaxies (e.g., Sc-Sd and Irr galaxies) and the fraction 
of mergers that host large initial MBHs becomes sub- 
stantially smaller at high redshifts. Future observations 
which would specifically search for high redshift dAGNs 
could test this prediction. 



In addition, we note here that the estimates of the 
dAGN LF and the separation distribution may be af- 
fected by some simplified assumptions made above. Be- 
low we discuss how our results would be affected by some 
other mechanisms, which may trigger nuclear activities, 
but is ignored in the model above. 

First, we have explicitly assumed that significant nu- 
clear activities can only be triggered in the gas-rich com- 
ponent of a major merger. If we relax this assumption 
and allow significant nuclear activities to be triggered 
in the gas-poor component of a major merger, then the 
estimated dAGN frequency would rise to about ten per- 
cent or higher at luminosity L[om] ^ 10 8 ' 5 £©, but is not 
affected much at lower luminosities. Under this assump- 
tion, the dAGNs at the high-luminosity end are primarily 
formed by major mergers of two red progenitor galax- 
ies as the initial MBHs in red galaxies are substantially 
big compared to those in blue galaxies with similar to- 
tal stellar masses. Mixed mergers of a red progenitor 
galaxy with a blue progenitor galaxy contribute little to 
the population of dAGNs with comparable luminosities, 
as the masses of the initial MBHs in these two progenitor 
galaxies are less likely to be comparable. 

Second, a dAGN can also be produced at larger sep- 
arations if nuclear activities could be triggered by tidal 
interactions when the two progenitor galaxies of a major 
merger are still far away from each other. Under this 
assumption, the model described in Section [2] can still 
be used but a large value of K should be chosen (e.g., 
K ~ 5 — 10). By assuming K ~ 5 — 10 and a con- 
stant Eddington ratio A = 0.25 till the completion of 
the mergers, our estimate of the frequency of kpc-scale 
dAGNs increases to ~ 10%. This increase is clearly be- 
cause most of the MBHs in the progenitor galaxies could 
now grow up significantly before the separation of the two 
components of a dAGN decays to the kpc scale. How- 
ever, increasing the value of the K parameter to ~ 10 
also implied that the predicted number of AGN pairs 
on scale of 10 — 100 kpc from our model is comparable to 
the total number of AGNs at luminosity L[oiii] ^ 10 7 -^©: 
which is not consistent with current observations. One 
possible way to solve this inconsistency would be to im- 
pose that tidal interactions induce only low-level nuclear 
activities (e.g., with A ~ 0.01 — 0.001) when the two pro- 
genitor galaxies are still far away from each other and 
the accretion rate is significantly enhanced only when 
the two progenitor galaxies become close enough to each 
other. In this case, a large number of paired AGNs on 
10 kpc — 100 kpc (with a frequency of several percent to 
ten percent) can be produced but the resulted frequency 
of dAGNs on kpc-scale is still the same as the current 
observations suggest. 

Third, successive mergers, in principle, can include 
mergers of two progenitor galaxies whose nuclei have al- 
ready been activated, which would also yield a dAGN. 
This scenario has not been discussed here in details 
mainly for simplicity. However, according to the ma- 
jor merger rates of galaxies presented in Section [21 the 
rate of successive major mergers, within the lifetime of 
detectable nuclear activities (e.g., < 1 Gyr), is extremely 
low, especially if successive merger events were assumed 
to be independent of each other. 

Future observations which would determine the distri- 
bution of dAGNs as a function of luminosity and separa- 
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tion may shed light on our understanding of the trigger- 
ing mechanisms, major mergers, tidal interactions and 
the mergers of galaxies with already-activated nuclei. 

4. CONCLUSIONS 

In this paper, we have constructed a phenomcnologi- 
cal model to estimate the number density of dAGNs and 
its evolution according to the observationally-estimated 
merger rates of galaxies and various scaling relations 
on the properties of galaxies and massive black holes. 
We have demonstrated that the observed frequency and 
separation distribution of dAGNs are compatible with 
the hypothesis that major mergers of galaxies lead to 
AGN/QSO activities and are consistent with most es- 
timates of galaxy major merger rates. We have also 
predicted that the frequency of those dual AGNs with 
double-peaked narrow line features becomes even smaller 
(i.e., 0.02% -0.06%) at higher rcdshifts z ~ 0.5-1.2 and 



the frequency of all kpc-scale dAGNs is - 0.04% - 0.1% 
at z ~ 0.5 — 1.2, which can be tested by future obser- 
vations. Future observations shall provide us with much 
better statistics on the distribution of kpc-scale dAGNs 
(and pairs of AGNs on even larger scale) as a function 
of luminosity and separation and hence will tighten the 
constraints on triggering mechanisms of nuclear activi- 
ties. 
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